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Triplet superconductivity in Sr2Ru04 is investigated with main interest on its internal 
degree of freedom. We perform a microscopic calculation to investigate how the chiral 
state d(k) = (k K ±ik y )z is realized among the underlying six degenerate states. Starting 
from the three band Hubbard model with spin-orbit interaction, we use a perturbation 
theory in order to calculate the pairing interaction. The p-wave superconductivity with 
T c ~ 1.5K is obtained in the moderately weak coupling region. It is shown that the or- 
bital dependent superconductivity (ODS) robustly appears in Sr2Ru04. We determine 
the stabilized state by solving the Eliashberg equations. We find that the Hund cou- 
pling term as well as the spin-orbit interaction is necessary for the "symmetry breaking 
interaction" . The main result is that the chiral state is stabilized in case of the p-wave 
symmetry with the main 7-band, which is obtained in the perturbation theory. When 
we assume the other pairing symmetry including the /-wave state, the symmetry break- 
ing interaction gives the other ci-vector. The electronic structure constructed from the 
^2 g -orbitals is essential for this result. 
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1. Introduction 

The discovery 1 ^ and establishment 2 ) of spin triplet superconductivity in Sr2Ru04 have stim- 
ulated much interests. Sr2Ru04 has the same crystal structure as high-T c cuprates and possesses 
a quasi-two-dimensional nature. 3 ' 4 ) In contrast to high-T c cuprates, 5,6 ) this system behaves as a 
typical two-dimensional Fermi liquid in the normal state. 3 ' 4 ) Three cylindrical Fermi surfaces are 
clearly observed by the quantum oscillation measurement. 3 ) These conduction band originates from 
the t2g-orbitals in Ru ions. 7 ' 8 ) The T 2 -law of the resistivity 4 ) and 3 ~ 4 times mass enhancement 3 ) 
indicate the importance of the electron correlation. The correlation effect is also suggested from 
the Mott insulating state in the related compound Ca2Ru04. 9 ) 

Theoretical interests have been focused on the unconventional superconductivity expected from 
the correlation effect and the low-dimensionality. The possibility of spin triplet superconductivity 
was first pointed out by Rice and Sigrist, 10 ) which was succeedingly supported by many experi- 
mental results. For example, NMR and NQR studies have shown no Hebel-Slichter peak, 11 ) the 
transition temperature (T c ) is suppressed by non-magnetic impurities, 12 ) and the time-reversal 
symmetry breaking is observed by /iSR. 13 ) The most significant evidence for spin triplet supercon- 
ductivity is obtained by the NMR Knight shift 14 ) and by the neutron scattering measurements, 15 ) 
* E-mail: yanase@hosi.phys.s.u-tokyo.ac.jp 
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which show that the spin susceptibility does not change through T c . The singlet pairing is denied 
by this fact. 

The triplet superconductivity has already been observed in heavy Fermion compounds, such as 
UPt3. 16 ) Probably, the recently discovered compounds UGe2 17 - ) and URhGe 18 ) are also the triplet 
superconductors. However, the microscopic investigation in these materials is generally difficult 
because of their complicated electronic structure. Actually, any microscopic theory acceptable for 
the superconductivity in heavy Fermion systems has not been constructed. On the other hand, the 
electronic structure of Sr2RuC>4 is relatively simple, that is, it has two-dimensional Fermi surfaces, 
a few degenerate orbitals, weak spin-orbit coupling and weak electron correlation. Therefore, 
Sr2Ru04 is the most precious and favorable compound for the microscopic investigation on the 
triplet superconductivity. Such studies on Sr2Ru04 will give various information also for the 
unconventional superconductivity in other materials. 

There have been several proposals for the microscopic calculation for the triplet supercon- 
ductivity in Sr2RuC>4, including the random phase approximation (RPA) 19 ' 20 ) and perturbative 
method. 21 ' 22 ^ In this paper, we use the perturbative method which will be justified in the moder- 
ately weak coupling region and be complementary to the RPA theories. If we use the spin isotropic 
(SU(2) symmetric) model, the possible triplet Cooper pairings are degenerate. In this paper, we 
take account of the spin-orbit interaction in the three band Hubbard-type model to lift the de- 
generacy, and show that the chiral state, which is consistent with experiments, is stabilized under 
the reasonable parameters (§3). In this case, the 7-band is mainly superconducting, and the k- 
dependence of the order parameter is classified into the p-wave symmetry (§2). It is shown that 
this pairing symmetry is essential to stabilize the d- vector along the z-axis (§4). 

The results obtained in this paper are related to the several interesting issues for S^RuO^; (a) 
direction of the ci-vector (b) pairing symmetry in fc-space (c) the power-law behaviors in various 
quantities, 23 ~ 27 ) and (d) mechanism of the triplet superconductivity. We briefly review these issues 
and the theoretical proposals for them. 

(a) d-vector. 

In the triplet superconductivity, the internal degree of freedom is an attractive character, i.e., 
the order parameter has three components in the spin space and is described by the <i-vector as 
A(fc) = \d(k)aa y . 28 ' 29 ^ However, the identification of the ci-vector using the microscopic Hamilto- 
nian is usually difficult, and actually not performed. Instead, it has been discussed phenomenolog- 
ically. 28 - 32 ) 

In case of Sr2Ru04, there exists a six-fold degeneracy together with the momentum component 
when the spin-orbit interaction is neglected. 32 ) The experimental results have concluded that the 
chiral state d(k) = (k x ± ik y )z is realized among the six-fold degeneracy, in which the time-reversal 
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symmetry is broken. 13 ' 14 ) * 

In order to explain why the <i-vector is parallel to z-axis, phenomenological calculations have 
been proposed by assuming the momentum dependence of the quasi-particle's interaction 33 ) or 
the anisotropy of the spin fluctuation. 34 ~ 36 ) In this paper, we carry out a microscopic calculation 
starting from the Hubbard-type Hamiltonian to identify the d-vector. The microscopic origin and 
the general properties of the "symmetry breaking interaction" are investigated. We find that the 
dominant contribution comes from the terms which have been overlooked in the phenomenological 
theories. 

(b) pairing symmetry. 

The issue on the pairing symmetry in fc-space is closely related to (c) the power-law behav- 
iors. 23 ~ 27 ) Because a finite excitation gap is expected in a pure p-wave chiral state, the under- 
standing of the gap-less behaviors is one of the main issues at the present stage. Actually, many 
theoretical studies have been dedicated to the resolution of this inconsistency. 

In order to compare various proposals, we explicitly write the momentum dependence of the 
order parameter. The chiral state is described as d(k) = (4> x {k) ± kj!>y(fe))i: where 4> x {k) (0 y (fc)) is 
a wave function having odd symmetry with respect to the reflection fc x — > — fc x (fc y — > — fe y ). For 
example, the tiny gap model has been proposed by assuming the p-wave function (j) x (k) = sinfe x . 37 ) 
Another candidate is the 2D /-wave symmetry; the f x 2_ y 2-wme state 38 ) (4> x (k) = sinfe x (cosfe x — 
cos fey)) and f xy -wme state 39 ) (4> x (k) = sinfe x sinfe 2 ) have been proposed. In relation with the 
thermal conductivity measurements, 40 ' 41 ) the 3D /-wave state with a horizontal line node has been 
proposed. 42-44 ) 

Although these proposals have succeeded in explaining the gap-less behaviors at least phe- 
nomenologically, a consensus remains to be obtained. In this paper, we compare the above gap 
functions including the orbital dependence by investigating whether chiral state is favored or not. It 
will be shown that the symmetry-breaking interaction given by the perturbation theory favors the 
chiral state only when the 7-band is mainly superconducting with p-wave symmetry. For example, 
another d- vector d{k) = 4> x (k)y ± 4> y (k)x is obtained in case of the f x 2_ y 2-wave symmetry. 

Recently, it was proposed that the power-law behaviors can be explained by the multi-band 
effect in a particular parameter set. 45 ) Then, three ideas play an essential role: (i) anisotropic gap 
structure in 7-band, 37 ) (ii) orbital dependent superconductivity (ODS) 46 ) and most importantly 
(hi) node-like structure in /3-band. 36 ' 45 ) Here, we use "ODS" as a situation where the amplitude 
of the order parameter strongly depends on the orbital. This proposal is consistent with the above 
argument on the ti-vector. Although we do not discuss this possibility, we find that the ODS 
naturally appears in Sr2Ru04. 

(d) mechanism of the triplet superconductivity. 

* Strictly speaking, this result relies on the /xSR experiment because the NMR Knight shift is measured under the 
parallel magnetic field. The NMR under the perpendicular magnetic field is desired for a stronger evidence. 
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The ferromagnetic spin fluctuations were first speculated 10 ' 47 ~ 49 ) in analogy with the superfluid 
3jj e 29,50) rp^jg p r0 posal, however, has been denied by the inelastic neutron scattering showing 
no sizable ferromagnetic spin fluctuations. 51 ) On the contrary, the neutron scattering detected 
the incommensurate spin fluctuations around q ~ (|7r, |7r), which were predicted in the band 
calculation 52 ) and in the multi-band RPA. 53 ) The incommensurate fluctuations have an anisotropic 
nature 54 ^ t Xz(q) > X±(q) which is derived from the spin-orbit interaction. 56 ) It has been proposed 
that this anisotropic spin fluctuations can produce the p-wave superconductivity with d-vector 
parallel to z-axis. 34 ~ 36 ) Since the incommensurate fluctuations are derived from the nesting of the 
a and (3 Fermi surface, this mechanism is favorable in the a- and /3-bands. 20 ) However, the more 
detailed calculation will be necessary since the 7-band has the largest density of states (DOS). 

Another proposal is the orbital fluctuations with the same wave vector q ~ (|7r, §7r). Takimoto 
has shown that the /-wave superconductivity is induced mainly on the a- and /3-band. 19 ) This 
mechanism requires that the inter-orbit repulsion is larger than the intra-orbit one. In general, this 
condition is difficult to be satisfied. 

Nomura and Yamada have proposed an another approach based on the perturbation the- 
ory. 21 ' 22 ) This approach is based on the general concept for the electronic mechanism. 57 ) The 
calculation for the single band model within the third order has given the p-wave superconductiv- 
ity. 21 ) It has been shown that the consistent results are obtained for the three band model where 
the 7-band is mainly superconducting. 22 ) 

In this paper, we similarly use the perturbation theory and show that the p-wave supercon- 
ductivity with T c ~ 1.5K is obtained in the moderately weak coupling region (§2.4). In this region, 
the third order term is sufficiently smaller than the second order term, but enhances the triplet 
superconductivity in cooperation with the second order term. The perturbative treatment for the 
"symmetry breaking interaction" consistently gives the chiral state. Thus, the triplet supercon- 
ductivity with appropriate (i-vector is obtained in the perturbative region when the characteristic 
band structure is correctly taken into account. 

Hereafter, we restrict the discussion to the unitary state in the 2D system. The fe z -dependence 
of the order parameter, and therefore the 3D f-wave state are not considered. Generally speaking, 
the 3D f-wave state 40-44 ) requires a strong fc z -dependence of the effective interaction. This case is 
not expected in the quasi-2D systems like Sr2Ru04. * 

The following part is constructed as follows. In §2, we explain the perturbation theory without 
spin-orbit interaction. The formulation is given in §2.1. The basic results including the ODS are 
explained in §2.2 and §2.3. In §2.4, the mechanism of the superconductivity and the justification of 
the perturbation theory are discussed. The effect of the spin-orbit interaction is taken into account 

t A query on this experimental result has been raised by the neutron scattering experiment. 55 ' 

* We estimate the three-dimensionality enough to give rise to the 3D /-wave state. The unrealistic value tx/ty ~ 0.1 
is required in the anisotropic cubic lattice. 58 - ) 
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in §3. The formulation is explained in §3.1. The results of the perturbation theory are shown in 
§3.2. The chiral state is stabilized when the 7-band is mainly superconducting. The comparison 
with the other paring states is carried out in §4. §5 is dedicated to the discussions on the validity 
of the approximations. The conclusion and some discussions are given in §6. 



2. Pairing Theory for the Three Band Model 

In this section, we discuss the spin triplet superconductivity in the multi-band Hubbard model 
without spin-orbit interaction. The perturbation theory is used to calculate the pairing interac- 
tion. 21 ' 22 ' 57 ' 59 ) The results are qualitatively consistent with Refs. 21 and 22. We explain the 
characteristic properties in the weak coupling region in more details. The robustness of the ODS 
is shown. The validity of the perturbation theory and the comparison with the other calculations 
are discussed. 



2. 1 Formulation 

The three band model adopted for Sr2Ru04 is described as, 

H = H + H LS + H h (1) 

3 

H = X^ £a ( fc ) C la, S C AVM> ( 2 ) 

a=1 k,s 

H LS = 2\Y,LiS t , (3) 

i 

H\ = U^^n ija ^n i ^ l + U'^ / ^n ij ani jb + JH^^(2S ij0 -S^ 6 + -n ij0 ra ij6 ) 

i a i a>b i a>b 

+ J Y1 Yl C i-k>,aA',a,f>°M c <l-kM> ( 4 ) 

a ^ b k,k',q 

= U /^2/^2 C q-fc',a,| C fc',a,T Cfc ' a 'T C 9-fc:aa 

a k,k',q 

+ ^ X/ C q-k',a,s C k',b,s' Ck ^ b > s ' C( l- k > a y' i 

a>b s,s> k,k',q 

+ ^ H XT5Z C q-k' ,a,s C l:' ,b,s' Cfe A sC< Z- fe > a ' s ' 

a>b s ,s' k,k',q 

+ J Y1 2 C l-k',aA',a,fk,b,^q-k,b,l- ( 5 ) 

a ^ b k,k',q 

The three bands are constructed from the three i2g-orbitals in the Ru ions. The index a represents 
the Ady Z (a = 1), 4<i xz (a = 2) and 4d xy (a = 3) orbitals, respectively. The kinetic energy term Hq 
is expressed by the two-dimensional tight-binding model, ei(fc) = — 2t' xy cosk x — 2t xy cosk y — fi xy , 
£2(fc) = — 2t xy cos k x — 2t' xy cos k y — fi xy and ss(k) = — 2i z (cos fe x + cos k y ) — 4t' z cos k x cos k y — /j, z 
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with t xy 3> t' xy and t z > t' z . Fig. 1 shows the typical Fermi surfaces which are consistent with 
the quantum oscillation measurement. 3 ) Hereafter, we describe the three bands without the spin- 
orbit interaction as x-, y- and z-band, respectively. The x- and y-bands have a one-dimensional 
nature and construct the a- and /3-bands through a weak hybridization. The weak hybridization is 
neglected for simplicity. The z-band has a two-dimensional nature and corresponds to the 7-band. 
Hereafter, we use the unit h = c = &b = 1 and choose the unit of energy as t z = 1. 




k x 

Fig. 1. Fermi surfaces with (solid lines) and without (dashed lines) the spin-orbit coupling. Parameters are chosen 
as t' z = 0.4, i xy = 1.5, = 0.2 and n x — n y — n z — 1.33. The spin-orbit coupling is chosen as A = 0.1. 

The spin-orbit interaction in the Ru ions are expressed by i?LS; where 2 A is the coupling 
constant. Typical value is considered to be A ~ 0.1 which is much smaller than the band width. In 
this section, we consider the case A = 0. 

The interaction term H\ represents the on-site Coulomb interactions including the intra-band 
Coulomb term U, inter-band Coulomb term U', Hund coupling term Jh and pair hopping term 
J. The parameters satisfy the relation U > 0, U' > 0, Jh < 0, and U > U' > \ Jh| ~ | J\ in the 
ordinary situation. The parameter J has an ambiguity in its definition which is determined by the 
relative phase between the orbitals. The standard notation gives a real positive value of J. 

Generally speaking, unconventional superconductivity originates from the momentum depen- 
dence of the effective interaction ^ Yl ^(k, k', a, a', b, b', s, s') c ! a s c -k a' s lC —k',V,s' c k',b,s- Here, we 
have only to consider the diagonal part V a; b(k, k' , s, s') = F(k,k',a,a,b,b,s,s') since the logarith- 
mic singularity appears from it. § Because the first order terms in the perturbation series do not 
contribute to the anisotropic superconductivity, the lowest order term is in the second order. In the 
following, we calculate all the second order terms with respect to Hi and the third order terms with 
the coefficient U 3 . This approximation is justified in the perturbative region U', \ Jh|, \J\ < U < W, 
where W is the band width W = 8 in our unit. We will explain the reason for including the third 
order terms in S2.4. 



We safely neglect the Cooper pairing between the different band. 
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In case of A = 0, the SU(2) symmetry in the spin space is conserved. Then, the effective 
interaction in the triplet channel is given by V* b (k,k') = V a ,b(k,k',s,s')\ s=s i regardless of the 
d- vector. Our procedure gives 

KbW) = V^(k,k') + VP(k,k')S a , b , (6) 
V™(k,k') = 2JJ HX ab(k-k') (for a + b), (7) 
V$(k,k') = -U 2 Xa(k-k')-(2U' 2 + 2U'J H + Jl)^2xa(k-k'), (8) 

a 

V^(k,k') = 2U 3 Re[^2G a (k + q )G a (k / + q)Ua(q)-UQm (9) 

q 

where G a (k) = (iuj n — e a (fc)) _1 is the Matsubara Green function. The indices k, k' and q are the 
four momentum k = (k, iu n ) etc., and the summation is defined as ^ fc = ^ Y2k ■ ^ ne index a in 
eq.(8) means the orbital with a ^ a. The functions Xa(q), Xab(l) and 4> a (q) are defined as, 

Xa(q) = -Y,Ga(k + q)G a {k), (10) 
k 

Xab(q) = -Re[J2G a (k + q)G b (k)}, (11) 
k 

Mq) = ^2G a (q-k)G a (k). (12) 

k 

The superconducting phase transition is determined by the linearized Dyson-Gorkov equation, 
A a (k) = -Y J Vlb(k,k')\Gb(k')\ 2 A b (k'), (13) 

b,k' 

where A a (k) is the anomalous self-energy on the a-band. The anomalous self-energy is an order 
parameter of the superconductivity. The momentum and orbital dependence of the anomalous self- 
energy represents the wave function of the Cooper pairs. The transition temperature is actually 
obtained by solving the Eliashberg equation, 

X c A a (k) = -^Vlb(k,k')\Gb(k')\ 2 A b (k'). (14) 

b,k' 

The maximum eigenvalue A e becomes unity (A c = 1) at the critical point. It is notable that the 
Eliashberg equation is reduced to the three independent equations if the off-diagonal term V^ b 
(oc JJh) vanishes. The transition temperature is uniquely determined owing to the pair hopping 
term J. 

Let us remark here the relation to the excitation gap below T c . The Eliashberg equation 
determines the momentum and orbital dependence of the anomalous self-energy at T = T c . Below 



T c , the Bogoliubov quasi-particle's energy is obtained as E a (k) = z a (k)^ e a (k) 2 + A£ x (fc) 2 where 
the renormalization factor z a (k) is given by z a (k)~ 1 = 1 — dReT^(k,(j)/du\ 0J= o and A° x (fc) = 
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\A R (k,E a (k))\, with A R (k) being the analytic continuation of A a (k). If we consider the weak 
coupling case, i.e. T c ,A^ x (fc) <C W, an approximation A R (fc, E a (k)) ~ A R (fc,0) ~ A a (A;,i7rT) is 
very precise around the Fermi surface. Therefore, the momentum and orbital dependence of the 
excitation gap is obtained from |A (fc,i7rT)| except for the factor arising from z a (k). Since we 
found that the a and fc-dependence of z a (k) is not outstanding in Sr2Ru04, the qualitative nature 
of the excitation spectrum below T c is obtained from the Eliashberg equation. Strictly speaking, 
the momentum and orbital dependence of A a (k) is deformed below T c . However, this deformation 
is usually small when T c <C W. 29 ^ 

2. 2 Order parameter and transition temperature 

In this section, the parameters and particle numbers are fixed to (t' z ,t xy ,t xy ) = (0.4, 1.5,0.2), 
(U,U f , Jh, J) = (5,1.5,-1,1) and n x = n y = n z = 1.33. The z-band has about 57 % of the 
total DOS under the above parameters. Using V£ b (k,k'), we obtain T c and the corresponding 
A a (k). In the summations, we divide the first Brillouin zone into 128 x 128 lattice points and take 
2048 Matsubara frequencies. We have confirmed that the obtained results do not depend on the 
numerical details, qualitatively. 

First of all, let us study the symmetry of the order parameter. The momentum and orbital 
dependence of A (fe) = A a (fc,i7rT) is shown in Fig. 2(a-c), where the temperature is chosen as 
T = T c = 0.0075. We find that the largest component is on the z-band (7-band). The pairing 
symmetry is the p-wave, which is consistent with Refs. 21 and 22. ^ The two nodes along k x = 
exist in the x and z-bands. There is no node in the y-band since the line k x = does not cross 
the Fermi surface. The mechanism of stabilizing the triplet superconductivity will be discussed in 
§2.4. 

From the symmetry, there is another solution 4> y (a,k) which is rotated 90 degrees from Fig. 
2. Generally speaking, arbitrary linear combinations A a (k) = pcfr x (a, k) + q(f> y (a, k) are degenerate 
at T = T c . We have shown the functions A a (k) = cf) x (a,k) in Fig. 2. The isotropic state like 
A a (&) = 4>x(a, k) ± i(/) y (a, k) is expected below T c in order to gain the condensation energy. 29 ) 

Second, let us see the angular dependence of the order parameter more closely. The amplitude 
of As(k) along the Fermi surface is shown in Fig. 3. Apparently A^(k) is small around (it, 0), which 
is qualitatively consistent with the proposal by Miyake and Narikiyo. 37 ) Although this behavior 
may be important for the gap-less behaviors observed experimentally, we should point out that the 
excitation gap in the chiral state is not so tiny as to explain the power-law behaviors. 

We would like to stress that the momentum and frequency dependence of A a (k) is determined 
so as to be most favorable for the superconductivity. In other words, the wave function of the 
Cooper pairs is spontaneously distorted so that the attractive interaction works most efficiently. 



In this paper, the pairing symmetry is defined by the sign reversal of the order parameter on the Fermi surface. 
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(a) x-band (b) y-band 




Fig. 2. Momentum dependence of the order parameter A a (fe) = A a (fe,i7rT). (a) a = 1 (cc-band), (b) a = 2 (jy-band) 
and (c) a = 3 (z-band). This figure shows the wave function x (a, fc). The signs + and — show the sign of the 
order parameter. The solid line in the basal plane shows the Fermi surface in each band. 
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Fig. 3. (a) Momentum dependence of the order parameter As(k) on the Fermi surface. The horizontal axis is the 
angle 8 = Arctan(fc y /fc x ) which is shown in (b). The results of the second order perturbation (long-dashed), and the 
third order perturbation at U — 3 (thick solid), (7 = 4 (dashed) and (7 = 5 (dash-dotted) are shown, respectively. 
Here, the inter-band interactions are chosen as U' = 0.3(7 and |Jh| = J = 0.2(7. The thin long-dashed line shows 
the result when only the third order terms are taken into account. The thin solid line shows the momentum 
dependence Aa(fc) = sinfc x . The wave function is normalized as the maximum value is unity (|A3(fc)| max = 1). 



It is interesting that the effective interaction V£ b (k,k') (eqs.(6-9)) is repulsive for the simplified 
wave function A a (k) oc sinfc x , but attractive for the distorted order parameter (shown in Fig. 2). 
Thus, the strong coupling theory based on the Eliashberg equation is essential for stabilizing the 
superconductivity. 
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Fig. 4. Transition temperature of the triplet superconductivity for t' z = 0.35 (triangle), t' z = 0.4 (circle) and t' z = 0.45 
(square). 

Third, we discuss the parameter dependences of T c . Figure 4 shows the n z - and ^-dependence. 
Apparently, T c increases with n z . This trend is consistent with the recent experiment. 60 ) We 
find that T c is not sensitive to the Van-Hove singularity. 61 ) The triplet superconductivity with 
T c ~ 0.007 is obtained around n z = 1.33 irrespective of the value t' z . If we assume W = 2eV, this 
critical temperature corresponds to T c ~ 20K. We have chosen the parameters where T c is higher 
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than the experimental value ~ 1.5K in order to keep the numerical accuracy. Actually, T c ~ 1.5K 
is obtained around U = 3.5 (see §2.4). We have confirmed that the qualitative features discussed 
in this paper do not depend on the values U and T c . 

2.3 Orbital dependent superconductivity 

It is important to note that the amplitude of the order parameter remarkably depends on 
the orbitals. The maximum value of the order parameter A a = |A a (/c)| max is considerably small 
for a = 1,2 (see Fig. 2). Thus, the superconductivity is dominated by the z-band. Such a strong 
orbital-dependence of the superconducting gap has been discussed phenomenologically 46 ) and called 
"orbital dependent superconductivity (ODS)". Hereafter, we denote the band in which the order 
parameter is largest as "main band" and the other bands as "minor bands" . The main band is the 
z-band in the present case. This is mainly because the DOS is largest in the z-band. 




Fig. 5. Ratio A3/A2 for U = 4 (circles) and {7 = 5 (triangles). 

We can see the robustness of the ODS behavior from Fig. 5 in which the i xy -dependence of 
the ratio A3/A2 is shown. As t xy decreases, the y-band becomes the main band since the DOS of 
the y-band increases. (Ai is always small if the symmetry is fixed to A a (k) = cf) x (a, k).) Figure 5 
shows that the decrease of A3/A2 is rapid, and the crossover region is very narrow. Thus, the order 
parameters in the different band only accidentally have the similar magnitude. The parameter t xy 
in the crossover region decreases with decreasing U. We find that z-band is the main band in the 
weak coupling limit, when t xy > 0.7. 

If the hybridization between the orbitals is strong, we expect the similar amplitudes of order 
parameters among the bands. In the present case, however, the coupling is only through the inter- 
band interaction V a ^ oc JJh which is much smaller than the intra-band one under the realistic 
condition. In other words, each band retains a character of the local orbital, and as a result, the 
ODS appears. This fact will be the basis of the calculation in §3 where we assume that the intra- 
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band interaction is most important. Even if the inter-band mixing between the x- and y-band is 
taken into account, the ODS behavior appears because the hybridization between the z-band and 
the other band vanishes in the 2D limit. This feature is robust owing to the inversion symmetry 
about the RuCVplane. 46 ) 

2.4 On the mechanism of the superconductivity 

From the present results for A = 0, we discuss the mechanism of the superconductivity. The 
following comment will justify the perturbation theory in this paper. 

Firstly, we wish to stress that the superconductivity is obtained within the second order per- 
turbation, i.e., in the lowest order perturbation. Furthermore, the results of the d- vector discussed 
in §3 do not change in the lowest order theory, qualitatively. Therefore, the obtained results are 
not the artificial ones in the third order theory. 

However, T c decreases exponentially in the small U region and thus we need to choose moderate 
values of U in order to obtain the experimental value T c ~ 1.5K. For example, T c ~ 1.5K is 
obtained around U = 3.5, U' = 0.3U and |Jh| = J = 0.2U, " where the expansion parameter 
for the perturbation theory is U /W ~ 0.38. We have confirmed that the absolute value of the 
third order term \V^ 3 \k, k')\ is still small in this region. For example, the average value of the 
third order term is about |W 3 )| ~ 0.4|W 2 )| in the above parameter set. (Here, we have defined 
the average value as \V^\ = f F f F \V$(k, k')\dkdk', where f F dk means the integration along the 
Fermi surface.) Therefore, we expect that the perturbation theory is qualitatively justified in this 
region. 

From the analysis below, we conclude that the role of the third order term is to stabilize 
the triplet superconductivity, i.e., significantly enhance the eigenvalue A e . This is because the 
momentum dependence is well suitable for the triplet superconductivity rather than the second 
order term. Therefore, the third order term plays an important role for the superconductivity in 
this region, which is the reason why we take into account the third order terms. 

To clarify the role of the third order terms, we compare the momentum dependence of As(fc) 
in Fig. 3 for various cases. If we consider only the second order terms, the order parameter on the 
z-band has a higher symmetry; it has 10 nodes on the Fermi surface. However, this behavior will 
be fictitious because, when the third order terms are introduced, the pairing symmetry becomes 
p-wave. This change is continuous, and the order parameter has similar form to the second order 
perturbation in the moderately weak coupling region. On the other hand, the qualitatively different 
form is obtained when only the third order terms are considered. These behaviors are natural 
because the system is in the perturbative region. 

Let us discuss here the results obtained in the perturbative calculations, RPA and fluctuation 

" Here, T c is determined by the extrapolation, because our numerical calculation is justified above T > 0.002 (~ 6K). 
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Fig. 6. The comparison of A e in various calculations for the triplet superconductivity. The solid line, circles, squares 
and triangles show the results of the second order perturbation (SOP), third order perturbation (TOP), RPA 
and FLEX, respectively. Here, the temperature is fixed to T = 0.005 (~ 15K) where A c ~ 0.4 corresponds to 
the experimental value of T c ~ 1.5K by the extrapolation. The parameters are chosen as U' = Jh = J — for 
simplicity. The dashed line shows the anti-ferromagnetic phase boundary U c = 2.0 in the mean field theory. 

exchange (FLEX) approximation. In Fig. 6, we compare the eigenvalue A e in the various calcula- 
tions for the triplet superconductivity. In the RPA, the observed T c ~ 1.5K is obtained only in the 
just vicinity of the magnetic order U = 2. This is mainly because the RPA term is not so suitable 
to the triplet pairing and the magnetic order is significantly overestimated in the RPA. It is also 
important that the contribution from the magnetic fluctuation to the triplet pairing is | of that 
to the singlet pairing. 49 ' 62 ) This fact can be remedied if the spin-orbit interaction is taken into 
account, 20 ' 34 ~ 36 ) but this issue is out of the scope of the present work. 

The magnetic order is actually suppressed by the fluctuations. In Fig. 6, we also show the 
results of the FLEX approximation 62,63 ) which is a convenient calculation to include the fluctua- 
tions. 64 ) We can see that the eigenvalue A e does not easily develop in the FLEX approximation. 
This is mainly because the de-pairing effect is very strong when the spin fluctuations have a local 
nature (namely, nearly qr-independent). ** Since the spin fluctuations arising from the z-band have 
a local nature, the spin fluctuation mechanism is generally difficult for the z-band. Actually, we 
find that the x- or y-band is the main band in the FLEX approximation for U > 2.5, where the 
pairing symmetry is the /-wave. 

These comparisons with the spin fluctuation theory indicate the importance of the non-RPA 
term which first appears in the third order correction. Actually, the third order terms significantly 
contribute to the enhancement of A e at U = 3, although their absolute value is still small. It is not 
easy to guess what happens in the higher order perturbation series. The actual result probably lies 
between the TOP and SOP results. In any case the triplet superconductivity is possible out of the 

** We should point out that the de-pairing effect is generally overestimated in the spin fluctuation theory. 59, 64 - ) This 
overestimation is serious when the spin fluctuations are spread in the momentum space. 
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on-site repulsive interactions, and the weak coupling treatment is meaningful. 
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Fig. 7. The comparison of A c for various symmetry. The symbols are the results of the TOP for triplet (circle), 
d x 2_ y 2-wave (square) and da^-wave superconductivity (triangle). The lines show the results of the SOP. The 
temperature is fixed to T = 0.005. The other parameters are chosen as U' = 0.3U and |Jh| = J = 0.2/7. 

Before closing this section, let us compare with the possible singlet superconductivity. Figure 
7 similarly shows the eigenvalue A e at T = 0.005. The results for the triplet, d x 2_ y 2-wscve and 
(i^-wave superconductivity are shown. The largest value of A c indicates that the corresponding 
state is stabilized in the lower temperature. If we consider only the second order perturbation, the 
d x 2_ y 2-wave symmetry is most stable. However, the third order term significantly enhances the 
triplet superconductivity, and slightly suppresses the d x 2_ y 2-wave superconductivity. We find that, 
although the RPA terms exist in the third order in the singlet pairing interaction V£ b (k, k'), they are 
considerably canceled by the non-RPA terms. This is an expected result when the spin fluctuation 
has a local nature. 64 ) As a result, the contribution from the third order terms to the d x 2_ y 2- 
wave pairing is small and negative. On the contrary, the third order non-RPA terms significantly 
enhances A e for the triplet pairing. Therefore, the triplet superconductivity is stabilized in the 
moderately weak coupling region (U > 1.5). 

The situation is very different for the case of High-T c cuprates. The cancellation in the third 
order is not strong in the d x 2_ y 2-w&ve channel because of the nesting property of the Fermi surface 
or the strong tendency to the anti-ferromagnetism. 64 ) Then, it is generally expected that the 
contribution from the non-RPA terms is not severe. Consequently, the spin fluctuation theory is 
qualitatively justified and the T c for the d x 2_ y 2 -wave superconductivity becomes very large. 

3. Determination of the D-vector 

In this section, we determine the d-vector assuming that the most effective spin-orbit interac- 
tions are those in the Ru ions. The SU(2) symmetry in the spin space is violated by this interaction, 
and thus the six-fold degeneracy of the triplet superconductivity is lifted. We clarify which pairing 
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state is stabilized by the spin-orbit interaction. 
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3. 1 Formulation 

First, we explain the formulation to take into account the spin-orbit interaction -£Tls- in this 
paper, H^s is included in the unperturbed Hamiltonian H' Q as 

H = H' + H h (15) 
H' = H + H LS . (16) 

Hereafter, we choose A = 0.1 unless we specify. The results do not depend on the choice as long 
as A is much smaller than the band width. By restricting the Hilbert space to the i2g-orbitals, the 
unperturbed Hamiltonian is written in the following way, 20 ' 33 ) 



k,s 





( £i(k) 


-sA 


- S A \ 


^ c k,l,s 
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c f c f c 1 " ) 


-sX 


e 2 {k) 
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Cfc,2,s 






^ -sA 


A 




\ Ck,3,-s 


) 



The sign s = 1 (s = — 1) corresponds to the up (down) spin. We denote the 3x3 matrix in eq.(17) as 
H' (k, s). Here, the wave functions of the t2g-orbitals have been chosen as 4d yz = i(4di +4d_i)/\/2, 
4d xz = — i(4di — 4d_i)/\/2 and 4d xy = — i(4c?2 — 4d_2)/v / 2, respectively. Note that a constant phase 
factor is multiplied to the 4d xz orbital in order to simplify the notation. This definition makes the 
pair hopping term J negative between the 4d xz and other orbitals. 

New quasi-particles are obtained by diagonalizing the matrix Hy(k,s) through the uni- 
tary transformation. This procedure corresponds to the transformation of the basis as 

( 4, M a U,s a k,3,s ) = ( c li,s c U,s c U,s )^( fc > s )> where C> ( fc ' s ) = Ki(M)) is a 
unitary matrix. New quasi-particles created by the operators a) kas are characterized by the pseudo- 
orbital a and pseudo-spin s. The pseudo-spin corresponds not to the spin but to the Kramers dou- 
blet. We can see that the rotational symmetry in the spin space is violated, and strictly speaking, 
the term "odd-parity superconductivity" should be used instead of the "spin triplet superconduc- 
tivity". Three Fermi surfaces of new quasi-particles have been shown in Fig. 1. The effect of the 
spin-orbit coupling on the band structure is weak, except for the mixing around k = (±|vr, ±|vr). 

The violation of the SU(2) symmetry typically appears in the spin correlation. We show 
the calculated results of the static spin susceptibility in the z-direction Xz(q) (Fig. 8(a)), and 
its difference from the planar component Xz{q) ~ X±{q) (Fig- 8(b)). We have used the Kubo 
formula and neglected the interaction term Hi, for simplicity. The static susceptibility has a 
peak around q = (|7r, |vr) which is consistent with the neutron scattering measurement. 51 ) The 
anisotropy results from the spin-orbit interaction. The z-component Xz{o) is considerably larger 
around q = (|7r, |7r) and slightly smaller around q = (0,0). These properties are also consistent 
with the NMR measurements 54 ) and the previous calculation. 56 ) 
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v ' ^ (b) X~X ± 




Fig. 8. (a) Static spin susceptibility of z-component Xz(<?)- (b) Anisotropy of the spin susceptibility X^io) ~ X±(l)- 
Here, t' z = 0.4, t xy = 1.25, t* y = 0.1 and A = 0.1. 

We have to calculate the effective interactions between the new quasi-particles. However, the 
interaction term Hi becomes very complicated by the unitary transformation. Therefore, we use 
the old basis Ck, a ,s tentatively, and at the end we apply the unitary transformation to obtain the 
effective interactions in the new basis. The effective interaction V between the old quasi-particles 
is calculated by the perturbative method within the same order as in §2. It should be noticed 
that the many terms are added to V through the off-diagonal Green functions. Moreover, we have 
to calculate the off-diagonal part of V because it contributes through the unitary transformation. 
Finally, the effective interaction T between the new quasi-particles is obtained as, 

T(k, k', a, a', /?', (3, si, s 2 , s 3 , s 4 ) = 
Y u *aA h ' s i) u *a>,a>(- k , s 2 )T(k, k' , a, a', b' , b, si,s 2 , S3, s 4 )u bl ^{-k', s 3 )u bi p(k', s 4 ), (18) 

a,a' ,b,b' 

and therefore, 

a,a',b,b' si,s 2 ,S3,S4 

= \ Y Y f( fe ' fc ^«'«^/ 3 ^/ 3 ' s l' s 2,S3,S4)a^ Sl a t _ fcIQ , )S2 a_ fc /^ IS 3a fc / I/ 3 IS4 . (19) 

a,a',/3,f3' s 1 ,s 2 ,s 3 ,s 4 

Here, we have redefined the up (down) spin in z-orbital as s = — 1 (s = 1). The diagonal part 
about the pseudo-orbital V a ^(k, k' , s±, s 2 , S3, S4) = T(k, k' , a, a, (3, (3, s\ , s 2 , S3, S4) is regarded as the 
pairing interaction. As a result, the Eliashberg equation is extended in the following way, 

\ c A a , Sl)S2 (k) = - £ V a: p(k,k\s 2 ,s 1 ,s 3 ,s i )\Gp(k')\ 2 A f3>S3:S4 (k'). (20) 

(3,k',S3,S4 
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Here, G a (k) = (ioJ n — E a (k))~ l is the Green function characterized by the pseudo-orbital where 
E a {k) is the energy of new quasi-particles. 

The practical calculation for the above equation is generally tedious since the pairing interaction 
V aj p(k, k' , si, s 2 , S3, S4) has so many terms. However, we can simplify the calculation by using the 
additional approximation based on the ODS argument combined with the perturbation with respect 
to the spin-orbit coupling A. As is shown in §2.3, the amplitude of the order parameter strongly 
depends on the orbital. In this case, the eigenvalue A e is almost determined by the interaction 
between the main orbital. This feature is justified for the pseudo-orbital since the spin-orbit 
coupling is weak. Therefore, it is sufficient to take into account only the diagonal part of the 
interaction V a (k, k', si, s 2 , S3, S4) = V a>a (k, k', si, S2, S3, S4), and investigate each case where the 
main pseudo-orbital is a = 1, 2 or 3. The validity of this approximation will be discussed in §5. 

The perturbation with respect to A is performed in the following way. First, the Green function 
G(k, s) = (io> n l — H' Q {k, s)) _1 and the unitary matrix U(k, s) are expressed as, 

( g 2 (k)g 3 (k) - A 2 -sX(g 3 (k) + A) -sX(g 2 (k) + A) \ 

-sX(g 3 (k)+X) gi (k)g 3 (k) - A 2 \(gi(k)+\) 
V -sA( 52 (fe) + A) X(gi(k) + \) gi{k)g 2 {k) - X 2 ) 
gi(k)g 2 (k)g 3 (k) - X 2 (g 1 (k)+g 2 (k)+g 3 (k)) - 2A 3 , 



G{k,s) 
A{k) 
U(k,s) 



/A(k), 



I 



\ 



z x (k) 1/2 

A 



A 



A 



e x (k)-£ y (k) 
A 



e y (k)-e x (k) 



\ 



e z (k)~e x (k) 
A 

e z (k)-e y {k) 



(21) 
(22) 
(23) 



e x (k)-£ z (k) e y (k)-£ z (k) 

where g a (k) = ito n — e a (k). The diagonal component z a {k) 1 / 2 is a normalization factor. Equation 
(23) is a correct expression within the lowest order about A, which is sufficient for the following 
calculation. It is important that the off-diagonal component of the matrices (eqs.(21) and (23)) 
appears in the order of A. From this observation, we find that the lowest order correction in 
the pairing interaction V a (k, k', s\, s 2 , S3, S4) is in the second order. Therefore, we perform the 
calculation up to 0(A 2 ) in the following. This procedure is justified because A is sufficiently small 
in Sr2Ru04 (X/W ~ 0.02). The normalization factor z a (k) = 1 — 0(A 2 ) also contributes within the 
second order. 

We find that the SU(2) symmetry in the d- vector space is violated by a partial set of the second 
order terms with respect to A, which we call the "symmetry breaking interaction". We find that 
the symmetry breaking interaction vanishes in case of Jh = 0. That is, the SU(2) symmetry in the 
d-vector space is violated by the combination of the Hund coupling and the spin-orbit interaction. 
Note that the SU(2) symmetry in the d-vector space can be conserved even if the SU(2) symmetry 
in the real spin is violated. 

As is shown in eq.(23), the expansion parameter is X/\e a (k) — £{>(&) |. This is sufficiently small 
(~ X/t z ) except for the crossing point of the Fermi surface around k ~ (±|7r,±|vr). We have 



18 Youichi Yanase and Masao Ogata 

confirmed that the contribution from this region is negligibie. We have also confirmed that the 
energy of the new quasi-particles, E a (k), is safely approximated as that of the corresponding old 
quasi-particles e a {k) which is the main component of the pseudo-orbital a. This approximation 
does not affect on the results within the numerical error. We define the a'-, f3'- and 7-band of the 
new quasi-particles which corresponds to the old x-, y- and z-band, respectively. Strictly speaking, 
a' and f3' are not the band, but constructed from the part of the a- and /3-band. Here, we have used 
the index representing the local orbital which is more convenient to discuss the superconductivity. 
The ODS means the strong dependence of the order parameter on the a'-, /?'- and 7-band. 

Finally, the eigenstates are classified by using the d- vector. From the symmetry argument, we 
have the following eigenstates (1) d(k) = k x x±k y y, (2) d(k) = k x y±k y x and (3) d(k) = (k x ±ik y )z. 
Here, k x and k y represent the reflection symmetry of the wave function and does not mean the 
detailed momentum dependence. There remains the two-fold degeneracy in each cases. Although 
other linear combinations are possible, we have chosen the symmetric states (l)-(3), which are 
expected to be stabilized in order to gain the condensation energy. The degeneracy in (1) and (2) 
is finally lifted by the weak mixing between the x- and y-orbitals, but the main results of this paper 
are not affected by this effect. A brief comment on this point will be given in §5. 

Because the symmetry d n (k) = 4> x (k) decouples with the different symmetry d n >(k) = 4> y {k) 
within our approximation, the effective interaction is defined according to the direction of the 
d- vector (n = x, y, z), 

V nj a(k, k') = ^(ia n a y )l l S2 V a {k, k' , s 1 ,s 2 , s 3 , S4){ia n a y ) S3:S4 . (24) 

Here, the spin part of the order parameter is assigned by the d-vector as, 

(*,,<*) A nW \ /-*(*) + 4(t) ) =iS(k) . (25) 

The resulting eigenvalue equation is obtained respectively, 

X e d n (k) = -J2 VnAk, k')\G a {k')\ 2 d n {k'). (26) 
k' 

It should be again stressed that the effective interaction V n ,a(k,k') depends on the main band a 
and the direction of the d- vector n. 

The numerical calculation is used to solve the three kinds of the Eliashberg equations (eq.(26)), 
which correspond to the three kinds of eigenstates (l)-(3). The three different transition temper- 
atures are obtained when Jh 7^ 0. We regard that the state with maximum T c is stabilized. This 
procedure is surely correct around T ~ T c . Strictly speaking, the determination of the ground 
state d- vector requires a direct estimation of the condensation energy, which was not performed in 
this paper. If the ground state d- vector is different from that for T ~ T c , a double transition is 
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indicated. However, there is no apparent reason for the double transition in the following sense. 
The condensation energy is maximumly gained in the above six eigenstates within the weak cou- 
pling theory for A = 0. We can show that the quadratic term in the Ginzburg-Landau expansion 
stabilizes the chiral state at finite A. Unless the higher order terms have an especial behavior, 
the stabilized state is retained below T c . Also we know that the double transition is not observed 
experimentally under the zero magnetic field. 23 ' 65 ) 

3. 2 Results of the perturbation theory 

First, we perform a fully microscopic calculation using the perturbation theory. In this case, 
the superconductivity is the p-wave and the main band is 7 as is shown in §2. This case is most 
feasible, since the DOS is largest in the 7-band. In this subsection, we show that the d-vector is 
also consistent with experiments. 

We calculate the diagonal part of the interaction between the 7-band, V n ^(k, k') and solve the 
corresponding Eliashberg equations (eq.(26)). Figure 9 shows the A-dependence of the transition 
temperature. We can see that T c is reduced by the spin-orbit coupling. This is mainly because the 
contribution from the zeroth order term Vs^k, k') monotonically decreases because of the normal- 
ization factor z z (k)z z (k') ~ 1 — 0(A 2 ). Although many terms arise from the spin-orbit interaction, 
they do not work as an attractive interaction. We confirm that T c at A = coincides with the 
result in the three band model (Fig. 4) within the numerical error. Thus, our approximation based 
on the ODS gives an appropriate estimation. 




Fig. 9. A-dependence of the transition temperature for n z = 1.33 and n z = 1.37. 

The temperature dependence of the eigenvalue A c for three kinds of the eigenstates is shown 
in Fig. 10. Hereafter, we choose the parameter set t xy = 1.25 and t' xy = 0.1 which increases the 
splitting of the degeneracy. We can see that the chiral state (3) d(k) = (k x ±ik y )z is stabilized in our 
calculation. We will show that this result is robust for the reasonable parameter set. Because the 
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splitting of the eigenvalue A e is almost independent of the temperature, we determine the stabilized 
state as a state with maximum eigenvalue A e at the fixed temperature T = 0.01, in the following. 

Note that the splitting of T c is considerably small, i.e., AT C ~ 0.01T C in this case. This 
is mainly because the symmetry breaking interaction is in the second order with respect to A, 
and requires the Hund coupling term Jh. The symmetry breaking interaction typically has the 
order ~ U JhA 2 /W 3 which is much smaller than the zeroth order term Vs^{k, k') ~ U 2 /W. Thus, 
the difference of the condensation energy is much smaller than the spin-orbit coupling energy A. 
Note that the spin-orbit coupling makes the spin correlation anisotropic even if Jh = (Fig. 8). 
Therefore, it is generally expected that the violation of the SU(2) symmetry in the <i-vector space 
is smaller than that in the real spin space. This result may be connected with the multiple phase 
diagram of UPt3. 30 ' 31 ' 66 ^ A brief discussion on this point will be given in §6. 




0.007 0.008 0.009 0.01 



0.007 0.0075 0.008 

T 

Fig. 10. Temperature dependence of A e for three kinds of the eigenstates (1) d(k) — k x x ± k y y (dashed- line), (2) 
d(k) = fc x y ± k y x (long dashed-line) and (3) d(k) — (fc x ± ik y )z (solid line). We can see that the state (3) is 
stabilized by the spin-orbit interaction. The inset shows the same quantities in the wide temperature region. 

The phase diagram is shown in Fig. 11, which is a main result of this paper. The chiral state 
(3) d(k) = (k x ± ik y )z is stabilized in the wide region around n z ~ 1.33. We have confirmed that 
the result in Fig. 11 is robust for the parameters of the a'- and /3'-band, only if the particle number 
is around n x = n y ~ 1.3. The other state (1) or (2) is stabilized only in the unphysical cases, for 
example, (i) hole-like Fermi surface of the 7-band (see Fig. 11), (ii) half-filled a'- and /3'-band, 
n x = n y ~ 1 or (iii) anti- ferromagnetic Hund coupling Jh > 0. Consequently, we conclude that the 
chiral state is robustly obtained in this calculation. 

4. Determination of the Pairing Symmetry 

In this section, we compare several pairing states from the view point of the internal degree of 
freedom. From the microscopic calculation in §3 we find that the symmetry breaking interaction 
and the main attractive interaction have a different origin. The superconductivity is mainly caused 
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Fig. 11. Phase diagram for the parameters of the 7-band. Circles and triangles represent the state (3) and (1), 
respectively. The state (2) is not stabilized. The solid line shows the number of particles n z = 1.33. The Fermi 
surface is electron-like (hole-like) in the left (right) side of the dashed line. The parameters for the a'- and /3'-band 
are chosen as (a) t xy = 1.25, t' xy — 0.1 and (b) f xy = 1.5, i xy = 0.2. 



by a part of the pairing interaction z a (k)V a (k, k')z a (k'), but this term does not contribute to 
the splitting of the degeneracy. Instead, the ci-vector is almost determined by the combination 
of the symmetry breaking interaction and the superconducting wave function </> x , y (&) at A = 0. 
This argument will be valid when the perturbation with respect to A is justified. In this section, 
we assume the pairing state at A = 0, and determine the (i-vector by calculating the symmetry 
breaking interaction using the formulation in §3. 

Before explaining the details, we summarize the obtained results in Table I where the stabilized 
state is shown for each pairing symmetry. We can see that the chiral state is stabilized only when 
the symmetry is the p-wave and the main band is 7. The other paring states give an inconsistent 
d-vector with the experimental result. Therefore, it is concluded that the most favorable pairing 
state is the p-wave on the 7-band, which is microscopically derived in the perturbation theory (§2). 



Table I. Obtained d-vector for each pairing symmetry and main band. 



4-1 P-wave symmetry on the ^-band 

First, we consider the pairing state obtained in §2 and §3. We assume that the total pairing 
interaction is given as 

V n „(k, k') = V<$(k, k') + V®(k, k'), (27) 

where vj£y(k, k') is the effective interaction derived from the microscopic Hamiltonian within the 
second order. We have neglected the third order terms (vjfy(k, k') oc U 3 ) because they do not con- 
tribute to the symmetry breaking interaction. Vy P \k,k') = z z (k)V^ (k, k')z z (k') is an attractive 
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interaction in the p-wave channel, where we choose V ( - p \k J k') = g(cos(k x — k' x ) + cos(/c y — k' y )). 
If g < — 1, the p-wave superconductivity is mainly caused by this term, and then we obtain the 
simple wave function 4> x {k) oc sinfc x . The symmetry breaking interaction is only from V^{k,k') 
and determines the d- vector. 
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Fig. 12. Phase diagram when the p-wave attractive interaction is assumed as eq.(27). Here, we choose g = —3. 
Triangles, squares and circles represent the state (1), (2) and (3), respectively. The solid and dashed lines are the 
same as in Fig. 11. 

The phase diagram is shown in Fig. 12 which is qualitatively the same as Fig. 11. Thus, the 
phenomenological model (eq.(27)) well reproduces the detailed calculation in §3. In particular, the 
chiral state (3) is stabilized around n z ~ 1.33. The only difference is that the state (2) appears 
in the region n z < 1.3 where the T c is very small in the perturbation theory We would like to 
stress that the chiral state (3) is robustly obtained independent of the detailed structure of the 
wave function 4> x (k). 



4-2 P-wave symmetry on the a- or (3-band 

Second, we discuss the d- vector when the symmetry is the p-wave but the main band is a or 
(3. This pairing state is also obtained by the perturbation theory when the DOS of a- and /3-band 
is large (§2.3). The notation a'- and /3'-band is convenient since the wave function (f> x (k) has large 
amplitude on one of them. When the symmetry is fixed to d n (k) = (f> y (k) (d n (k) = (f) x (k)), the main 
band is a' {(3'). Since the a'-band is related to the /3'-band through the 7r/2-rotation, d(k) = k y x 
on the a'-band is degenerate to d(k) = k x y on the /3'-band within our approximation. Therefore, 
the same classification (1), (2) and (3) is valid. 

Since this pairing state is microscopically obtained, we show the results of the perturbation 
theory which has been formulated in §3. We similarly calculate the diagonal part of the interaction 
Vn,(3'(k,k') (or V n ,a'(k,k')). Then, the wave function on the /3'-band has a similar momentum 
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dependence with Fig. 2(b). 

We find that the state (2) d(k) = k x y ± k y x is stabilized under the reasonable parameter 
set, which is inconsistent with the ^SR experiment. The splitting between the states (2) and (3) is 
smaller by an order of magnitude because the term with coefficient UJuX 2 vanishes in the symmetry 
breaking interaction in this case. Instead, the splitting energy has a coefficient U'JuX 2 or J^A 2 . 
Therefore, it is conceivable that the stabilized state changes depending on the details neglected 
in our approximation. However, the weak hybridization between the x- and y-orbital generally 
stabilizes the state (1) or (2) as will be discussed in §5. It is therefore expected that the chiral state 
(3) is not favored in this case. 

These results are robust in the wide parameter region, 1.0 < n z < 1.6, 1.26 < n x = n y < 1.64, 
and 0.05 < t' xy < 0.5. Even if the DOS of the a'- and /?'-band is very large (t xy = 0.7 and 
t' xy = 0.056), the stabilized d- vector does not change. The qualitatively same results are obtained 
when we assume the attractive interaction as in eq. (27) . 

We comment on another possibility of the p-wave superconductivity where the main band is a' 
(/?') for the wave function d n (k) = (f) x (k) (d n (k) = 4> y (k)). In this case, nodal quasi-particles exist 
in all of the states (l)-(3), and power-law behaviors 23-27 ) can be explained. We similarly investigate 
this possibility by assuming the pairing interaction like eq.(27) and conclude that the d- vector is not 
parallel to z under the reasonable parameter set (see Table. I). Therefore, this candidate (p-wave 
symmetry with node) is also inconsistent with the time reversal symmetry breaking. 13 ) 

4-3 F-wave symmetry 

Finally, we investigate the case of the /-wave symmetry which has gap-less quasi-particles in 
the chiral state. For this candidate, we assume again a phenomenological pairing interaction like 
eq.(27) and choose V^ p \k,k') = g(4> x (k)4> x (k') + <fr y (k)<j) y (k')). The / x 2_ y 2-wave state is obtained 
by choosing (p Xjy (k) = sin k X:Y (cos k x — cos/c y ). The eigenstates of the Eliashberg equation are 
similarly classified into (1') d(k) = (k x - k y )(k x x ± k y y), (2') d(k) = (k x - k y )(k x y ± k y x) and (3') 
d(k) = {k 2 x -k 2 y ){k x ±\k y )z. 

First, we consider the case where the main band is 7. In this case, the nodal quasi-particles 23-27 ) 
and the time-reversal symmetry breaking 13 ) coexists in the chiral state (3'). From the comparison 
with the thermal conductivity data, it has been proposed that the state (3') is the best candidate 
for S^RuOzt. 38 ) However, our result shows that the realized state is (2') for n z > 1.26 and (1') for 
n z < 1.26. The state (3') is not stabilized by the spin-orbit interaction. Thus, this candidate is not 
probable at least in the weak coupling region. 

The other candidate (3") d(k) = k x k y (k x ± ik y )z has been proposed by Graf and Balatsky. 39 ) 
However, the stabilized state is not (3") but (1") d(k) = k x k y (k x x± k y y) in our results. Thus, the 
chiral state is not stabilized also in the / xy -wave state. 
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Next, we discuss the case where the main band is a or [3. The candidate is d n (k) = 
sin fc y (cos k x — cos k y ) on the a'-band (or equivalently, d n (k) = sin fc x (cos k x — cosfc y ) on the [3 1 - 
band). We have found that this pairing state is most favorable in the triplet channel when t xy /t z 
and U/W is sufficiently small, or when the RPA or FLEX approximation is applied for small t xy /t z 
(§2). We have confirmed that the stabilized state is not (3'), but (f ') d(k) = (k x — k 2 )(k x x ± k y y). 
Therefore, this candidate is also incompatible with the /j,SR experiment. In Takimoto's pro- 
posal based on the orbital fluctuations, the main band is a' or /3', and the wave function is 
approximately written as A a r(k) = smk y (A + B(cosk x + cosk y )) (A = 0.51 and B = 1.3). 19 ) 
We have also carried out the calculation for this pairing state and confirmed that the state 
d(k) = (A + £>(cos k x + cos fc y ))(sin k x x ± sin k y y) is stabilized. 

The above results have been summarized in Table I. It should be understood that the chiral 
state is stabilized only in the particular case. 

5. Discussions on the approximations 

We have used the following approximations in order to determine the d- vector. 

(i) Perturbation on the electron correlation. 

(ii) Perturbation on the spin-orbit interaction. 

(iii) Contribution from the minor band has been neglected. 

(iv) Hybridization between the x-and y-orbitals has been neglected. 

In this section, we show that the above approximations are justified under the respective conditions 
which are expected for Sr2Ru04. 

First, the approximation (i) is a main assumption of this paper; the moderately weak interac- 
tion is assumed. Conventional Fermi liquid behaviors 3 ' 4 ) have supported the weak coupling theory. 
It is necessary for the argument in §4 that the perturbative treatment for the "symmetry break- 
ing interaction" is justified. This assumption is more promising because the symmetry breaking 
interaction requires the Hund coupling term, and is expected to be small. 

Second, the approximation (ii) is safely justified since the spin-orbit interaction is sufficiently 
weak in Sr2Ru04 (2A ~ O.leV) compared with the band width (~ 2eV). This approximation is 
probably not justified in the heavy fermion compounds. 

Third, we discuss the contribution from the minor band which is negligible in case of the 
ODS. 46 ) In the starting model at A = 0, the order parameter in the different band only weakly 
couples through the inter-band pairing interaction, V a ^(k,k'). Therefore, the order parameter in 
the minor band has much smaller amplitude than that in the main band. This situation is similar 
to the two-band superconductivity in the s-wave case. 67 ) The ratio furthermore decreases when 
the intra-band pairing interaction in the minor band is repulsive. The contribution from the minor 
band to the symmetry breaking interaction has at most an order |V^ 5 b|A 2 /iy 2 oc JJu\ 2 /W 3 , which 
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is smaller than the main contribution (~ UJftX 2 /W 3 ). Thus, the approximation based on the ODS 
is justified at finite A. 

Finally, we discuss the weak hybridization between the x- and y-orbitals. Taking the hybridiza- 
tion into account, the unperturbed Hamiltonian eq.(17) is changed as, 

^ £i(fe) ig(k)-s\ -sX \ ( c fe; i iS ^ 

-ig(k) - sX e 2 {k) X c k ^ s ■ (28) 

\ sX X e 3 (k) J \ c fci3 _ s J 

The tight binding model gives g(k) = At m sin k x sin ky , where t m is the hopping integral between 
4d xz - and 4<i yz -orbitals in the next nearest-neighbor site. We can treat this mixing term in a 
perturbative way as is used for the spin-orbit interaction. That is, the perturbation with respect 
to the off-diagonal term in eq.(28) is possible. Two kinds of the effective interactions newly appear 
within the second order perturbation, namely the terms proportional to (a) and (b) it m X. The 
term (a) is not important since it does not violate the SU(2) symmetry. We discuss the cross term 
(b) oc it m A which contributes to the symmetry breaking interaction. 

We easily find that the term (b) vanishes when the main band is 7, and therefore, the ap- 
proximation (iv) is well justified. On the contrary, the term (b) appears when the main band 
is a' or (3' . The pairing state with different reflection symmetry (for example, d(k) = k x x on 
the a'-band and d(k) = k y y on the /?'-band) couples through the term (b). Therefore, the 
term (b) mainly works through the inter-band interaction V^^k, k' , s±, S2, S3, S4) in the present 
case. This interaction has a finite value only when s± = S2 = S3 = S4, and satisfies the 
relation V^,(k, k',s) = V^,(k, k' , s, s, s, s) = -V^,(k,k' - s) and V^iK, k y , k' x , k' y , s) = 
— V^, b ^,(— /c x , k y , —k'y,, k' y , s) = — V^, b ^,(A;x, — k y , k' x , —k' y , s). Consequently, we find that the two-fold 
degeneracy in the eigenstates (1) and (2) are lifted by the term (b). The other effect is neg- 
ligible. Then, we obtain the eigenstates (la) d(k) = k x x + k y y, (lb) d{k) = k x x — k y y, (2a) 
d(k) = k x y + kyX and (2b) d(k) = k x y — k y x without the degeneracy. These eigenstates correspond 
to the one-dimensional representations in the symmetry argument. 28 ) On the other hand, the chiral 
state (3) is not affected by the term (b), and therefore two-fold degeneracy is maintained. We find 
that one of the states between (la) and (lb) ((2a) and (2b)) has a lower energy owing to the term 
(b), and the other has a higher energy. This is because the matrix element arising from the term 
(b) has an opposite sign for them. Consequently, two of the states (la), (lb), (2a) and (2b) are 
favored by the term (b), and accordingly the state (3) is unfavored. In other words, the d- vector 
along the z-axis is not stabilized by the term (b). Since the term (b) is in the first order with 
respect to A, the state (3) is not stabilized in the weak coupling limit about A when the main band 
is a or [3. 

As a result, the approximation (iv) does not affect the obtained results in Table I. The mixing 
term works only when the main band is a or (3. Then, the state (3) is not stabilized in any case. 
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6. Summary and Discussion 

In this paper, we have theoretically investigated the internal degree of freedom of the triplet 
superconductivity in Sr2Ru04. The obtained results are summarized in the following way. 

(1) The perturbation theory gives a p- wave superconductivity in the moderately weak coupling 
region. The ODS is a robust result when the inter-orbital interaction is smaller than the intra- 
orbital one. Then, the order parameter has an exceedingly large amplitude in one of the band a ! , 
P' or 7. The main band is 7 under the reasonable parameter set. 

(2) The spin-orbit interaction in Ru ions makes the <i-vector parallel to the z-axis. In this case, 
the chiral state d(k) = (k x ± ik y )z with time-reversal symmetry breaking is expected. The orbital 
and momentum dependence of the order parameter is essential for this result. 

(3) When the other pairing state is assumed, the (i-vector is not pointed to the i-direction. 
The /-wave pairing state and the pairing state with main a'- or /3'-band are included in such cases. 
These cases are not consistent with the time reversal symmetry breaking indicated by the /iSR 
experiment. 13 -* 

(4) The splitting of T c is small even if the spin-orbit interaction is much larger than the 
characteristic energy of the superconductivity, namely A 3> T c , A. If we put the parameters as 
2A = O.leV and W = 2eV, the splitting is about AT C ~ 0.04T C . This is mainly because the 
"symmetry breaking interaction" is in the second order with respect to X/W. The first order 
term vanishes because the band mixing does not exist at A = 0. It is an another reason that the 
symmetry breaking interaction requires the Hund coupling term. Since |Jh| is smaller than the 
intra-band repulsion U, the symmetry breaking interaction has a smaller coefficient than the main 
part of the effective interaction. 

From the results (l)-(3), the present work has supported the weak coupling theory for 
Sr 2 Ru0 4 . 21 ' 22 ) It has been confirmed that the perturbation theory gives the probable pairing 
state, namely the p-wave symmetry and the main 7-band, and furthermore the consistent d- vector 
along the z-axis. 

Then, the detailed analysis of the gap-less behaviors 23-27 ^ will be important for the coherent 
understanding. If the discussion is restricted to the single band model, the line node appears only 
accidentally. This situation is not improved even if we consider the three-dimensional model 42-44 -* 
because the line node is not necessarily horizontal. Only the point node is allowed from the 
symmetry argument. 68 ) We think that this difficulty can be resolved by the multi-band effect 
within the two-dimensional model. Some proposals have been given along this line. 36 ' 45 -* We 
consider that the low-energy excitation is mainly included in the a- and /3-band, while T c and 
(i-vector are determined by the 7-band. To show the gap-less behaviors, an detailed investigation 
on the minor bands is essential. Such study is beyond this paper and will be an important future 
issue. 
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Here, we would like to point out that our systematic treatment for the spin-orbit interaction 
is necessary to discuss the d- vector. The spin-orbit interaction determines the d- vector through 
(i) the virtual process in the effective interaction and (ii) the unitary transformation of the quasi- 
particles. When the chiral state is stabilized, the dominant contribution comes from the cross-term 
of (i) and (ii). Therefore, we have to treat the two effects (i) and (ii) in a same footing. The 
estimation for the effect (i) only 20 ) or the effect (ii) only 33 ) is inadequate for the microscopic 
theory. We should further mention that the effective interaction derived here is quite different from 
the phenomenological assumption in Ref. 33. Indeed, the chiral state is not stabilized if we take 
into account only the effect (ii). 

At last, we comment on the phase diagram of UPt3 from the result (4). A heavy Fermion 
system UPt3 is a spin triplet superconductor, and shows three different superconducting phases 
under the magnetic field. The phase diagram has been explained by assuming the weak spin-orbit 
interaction 30 ^ because this assumption is necessary to explain the NMR Knight shift. 66 ) However, 
this assumption has raised a serious question since the spin-orbit coupling is generally strong in the 
heavy Fermion system. 31 ) We think that the result (4) gives a clue to this question. It has been 
shown that the Hund coupling term is needed to violate the SU(2) symmetry in the d-vector space. 
On the contrary, the SU(2) symmetry in the real spin space is violated even if Jh = (see Fig. 8). 
Therefore, the violation of the SU(2) symmetry can be much smaller in the d-vector space than in 
the real spin space. 69 ) In other words, it is possible that the anisotropy is almost absorbed in the 
character of the quasi-particles, and only a weak anisotropy is remained in the residual interaction. 
Thus, the microscopic theory in this paper gives a general insight on the triplet superconductivity. 
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